Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



O 
O 

(N 

5-1 . 

< 

ON 



6 

CO 



THE DEPENDENCE OF THE ESTIMATED LUMINOSITIES OF ULX ON SPECTRAL MODELS 
A. Senorita Devi 1 , R. Misra 2 , V. K. Agrawal 3 and K. Y. Singh 1 

Draft version February 5, 2008 

ABSTRACT 

Data from Chandra observations of thirty nearby galaxies were analyzed and 365 X-ray point sources 
were chosen whose spectra were not contaminated by excessive diffuse emission and not affected by 
photon pile up. The spectra of these sources were fitted using two spectral models (an absorbed 
power-law and a disk blackbody) to ascertain the dependence of estimated parameters on the spectral 
model used. It was found that the cumulative luminosity function depends on the choice of the 
spectral model, especially for luminosities > 10 40 ergs/s. In accordance with previous results, a large 
number (~ 80) of the sources have luminosities > 10 39 ergs/s (Ultra-Luminous X-ray sources) with 
indistinguishable average spectral parameters (inner disk temperature ~ 1 keV and/or photon index 
r ~ 2) with those of the lower luminosities ones. After considering foreground stars and known 
background AGN, we identify four sources whose minimum luminosity exceed 10 40 ergs/s, and call 
them Extremely Luminous X-ray sources (ELX). The spectra of these sources are in general better 
represented by the disk black body model than the power-law one. These ELX can be grouped into 
two distinct spectral classes. Two of them have an inner disk temperature of < 0.5 keV and hence 
are called "supersoft" ELX, while the other two have temperatures > 1.3 keV and are called "hard" 
ELX. The estimated inner disk temperatures of the supersoft ELX are compatible with the hypothesis 
that they harbor intermediate size black holes, which are accreting at ~ 0.5 times their Eddington 
Luminosity. The radiative mechanism for hard ELX, seems to be Inverse Comptonization, which in 
contrast to standard black holes systems, is probably saturated. Extensive variability analysis of these 
ELX, will be able to distinguish whether these two spectral class represent different systems or they 
are spectral states of the same kind of source. 
Subject headings: Galaxies: general - X-rays: binaries 
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1. INTRODUCTION 

In the last few years, Chandra observations of nearby 
galaxies have detected many non-nuclear X-ray point 
sources jK aaret e t ali [2001 1 : iMatsumoto et all 120011 : 
IZezas fc Fa bbiano 2002( 1 , some of which have isotropic lu- 
minosities > 10 39 ergs/s and are called Ultra luminous X- 
ray Sources (ULX). While some of these ULX are super - 
nova remnants (e.g. iRvder et all 119931 : iFox et al.lf2000h . 
it is believed that the majority of them are compact ac- 
creting systems. Indeed, ASCA X-ray spectral studies of 
many ULX have revealed that t hey display the charac- 
teristics of accretin g black holes (Makish ima et al"]l2000L 
iMizuno et al.ll2001l ) . ULX have also bee n called Super- 
Eddington Sources ( Fabbiano 1989, 2004) and Intermedi- 
ate luminosity X-ray object s ((Roberts fc Warwick! 120001 : 
iColbert fc Mushotzkvlfl99"9h . 

Since these ULX sources emit radiation at a rate 
larger than the Eddington luminosity for a ten-solar mass 
black hole, they are believ ed to harbor a black hole of 
mass 10 M. <M < 10 5 M r , ((Colbert fc Mushotzky|[l999l : 
iMakishima et alj |2000|) where the upper limit is con- 
strained by the fact that a more massive black hole would 
have settled into th e nucleus due to dynamical friction 
(|Kaaret et qf J[200TT ). Black holes in this mass range are 
called Intermediate Mass Black Holes (IMBH), since they 
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seem to represent the missing component of the black 
hole mass spectrum with masses prevailing in the gap be- 
tween those of stellar mass black holes found in Galactic 
X-ray binaries and those associa ted with Active Galac - 
tic Nuclei, M ~ 1 s - 10 9 (iRichstone et al.| [l998h. 
iMiller fc Colbert! (|2004h and iMillerl (120051) review the 
presen t evidence for IMBH in ULX and iLiu fc Mirabel 
( 2005) have compiled a catalogue of some ULX and prop- 
erties. 

Alternate models for ULX are that th eir luminosi- 
ties are super-Eddington (|Begelmanl I2002T ) or that their 
emission i s beamed from a geometrically thick accre- 
tion disk ((King et a/.l (200D. However, it has been ar- 
gued that in the latter case, such thick "funnel" shaped 
disks enhance the obse rved flux by just a factor of few 
((Misra fc Sriraml |2003T ). For all of these models, the 
creat i on of such sources dPortegies. Zwart fc McMillian 
I2002t iTaniguchi etall l2000t iMadau fc Reesl l2001h and 
process by whic h they sustain high accretion rates 
(King et al. 2003), are largely unknown. 

Investigations on the nature of ULX have been un- 
dertaken by studying the spectra and variability of in- 
dividual sources. For e xample analysis o f the spectra 
of NGC 1313 X-l. X-2 (IMiller et aLl l2003h and M81 X- 
9 ((Miller. Fabian fc Millerl |2004|) . revealed the presence 
of a cool accretion disk component (kTi n ~ 0.1 — 0.5 
keV), suggesting that ULX indeed harbor IMBH. Tran- 
sitions between two spectral states, similar to those 
seen in Galactic b lack hole systems, have been reported 
in NGC 1313 X-l dColbert fc Mushotzkvl 1 19991 1 and two 
sources m IC342 ((Kubota et q/.ll2001h . Spectral transi- 
tions have also been reported in two sources in NGC 1313 
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(|Feng fc Kaaretll2006D . 

The large collecting area of XMM-Newton, allows for 
detailed spect ral fits to ULX, wh i ch often comprise of tw o 
compone nts (IWang et all 12004 iFeng fc Kaaretl l2005f ). 
However, iGoncalves fc Soria (2006), have argued that 
such soft spectral components depend on the complex- 
ity of the fitting model. An interesting object is the 
brightest X-ray point source in M 82, whose intrinsic 
luminosity has been measured to be as high as 1.6 x 
10 41 ergs/s ()Ptak and Griffiths! Il999l) . The detection of 
a 54 mHz quasi-periodic oscillation in its X-ray light 
curve suggests that the source is a compact object and 
not a background AGN (jStrohmaver fc Mushotzkvll200l 
iDewangan et al.ir2006al ). The spectra of this source can 
be fitted by a power-law with photon index, r ~ 2 
(|Fiorito fc Titarchuk|[200l . but is more consistently fit- 
ted with a flatter power-law with an high energy cut- 
off around ~ 6 keV, which can be in terpreted as opti- 
cally thick, saturated Comptonization (|Agrawal fc Misral 
2006). A quasi-periodic oscillation has also been discov- 
ered in the bright X-ray source of Holmberg IX, which 
is similar to the source in M82 in havi ng a flat spec- 
trum (r ~ 1) with a ~ 9 keV cu toff llDewangan et al.1 
2006b). Recently. IStobbart et all (|2006l) found that the 
XMM-Newton spectra of eleven of the eighteen ULX 
studied by them, showed such high energy curvature. 
Chandra observations of NGC 5204 X-l, also reveals the 
presence of an optic ally thick Comptonized component 
([Roberts et al.1 l2006f ) . While these results of individual 
ULX are intriguing, there does not seem to be any signif- 
icantly distinguishable spectral property of ULX, and in 
general their spectra can be described either by steep or 
flat power-law indices, with and without soft components 
(e.g. IDewangan et al.ll200l . 

Another line of investigation is to construct the cu- 
mulative luminosity function and histograms of spec- 
tral parameters of a l arge sample of X-ray sources (e.g. 
IColbert fc Pt ak 2002). The hope here is that, in case 
ULX are a distinct class of sources and/or they can be 
classified into distinct subgroups, the luminosity func- 
tion should exhibit a break a nd their spec t ral pa rame- 
ters should show clustering. ISwartz et all ((2004) ana- 
lyzed data from 82 galaxies and estimated the luminos- 
ity function and spectral parameters of the X-ray point 
sources in these galaxies. They found that the average 
photon index (as well as the distribution) of the ULX 
and the less luminous sources is nearly same. Moreover, 
their spatial and variability distributions are also similar. 
While, their analysis revealed that the luminosity func- 
tions of ULX depend on the host galaxy type and star 
formation rate, they did not find any significant evidence 
for breaks in them. 

As emphasized bv lSwartz et al\ (|2004h . the power-law 
model they used to fit the data was chosen as an empir- 
ical one. They attempted to fit all the sources with the 
power-law model and only for those sources that did not 
provide a reasonable fit, they used other models like disk 
black body. They note that for many sources that are 
well fitted by the power-law model, other spectral models 
can also represent the spectra. In this work, we consider 
a smaller sample of 30 galaxies but fit the spectra of the 
points sources, with both a power-law and a disk black- 
body model. In principle, the intrinsic (i.e. the absorp- 
tion corrected) luminosity inferred for a source may be 



different for the two spectral models. Our first motiva- 
tion here is to make a qualitative estimate of this differ- 
ence by noting the dependency of the luminosity function 
on the spectral model used. A second motivation for this 
work is based on the expectation that estimations of a 
different spectral parameter like the inner disk tempera- 
ture (as compared to the photon index), maybe better in 
distinguishing ULX from other sources or they may re- 
veal dependencies (like correlations between luminosity 
and temperature) which could shed light on the nature 
of these sources. As mentioned earlier, the spectra of 
several ULX are complex requiring more than one com- 
ponents. Thus, the two spectral models, an absorbed 
power-law and and an absorbed disk blackbody, should 
be considered as empirical ones, which can adequately 
represent low count data and these models need not be 
the correct physical representation of the actual source 
spectrum. 

It is clear that U LX are fairly commo n in nearby galax- 
ies. For example, ISwartz et al\ (|2004l ) identified 154 of 
them in 82 galaxies. If these sources are Eddington lim- 
ited, then black holes > 1OM may indeed be quite com- 
mon. However, only a few of these sources have luminosi- 
ties greater than 10 40 ergs/s. It is these sources, (the 
best example being the bright source in M82 X-l), that 
require black hole with masses, M ~ 10 2_4 M Q . The de- 
velopment of a self-consistent theory which explains the 
process by which such black holes are created and un- 
dergo high accretion, is theoretically challenging. Hence, 
it is important to estimate the number of sources whose 
minimum intrinsic luminosity exceeds 10 ergs/s. Here 
the minimum value of the luminosity should not only in- 
clude the statistical spectral fitting error, but also the 
variations in the luminosity estimation that may occur 
upon using different viable spectral models. The third 
motivation of this work is to identify such sources which 
we call Extremely Luminous X-ray sources (ELX). To 
avoid possible ambiguities in determining the luminos- 
ity, we have chosen only those sources whose spectra are 
not contaminated by excessive diffuse emission and which 
are not affected by photon count pile up. Identification 
of such relatively "clean" systems would allow for more 
detailed studies of their properties to be undertaken and 
would be the first step toward understanding their na- 
ture. 

A similar a nalysis has been undertaken recently by 
iWinter et al.l (|2006D for thirty galaxies observed by 
XMM-Newton. Since their motivation was to check if 
ULX also exhibit soft and hard spectral states like black 
hole binaries, they limited their analysis to bright sources 
where detailed spectral analysis could be done. They 
found 16 sources as possible low-state ULX and 26 as 
high state ones. For the high state ULX, the observed 
range for the black body temperature was 0.1 — 1.0 keV, 
with the more luminous sources having lower temper- 
ature. The Chandra analysis undertaken here, can be 
compared with this contemporary work and as we discuss 
later, the analysis are more or less consistent with each 
other. This consistency is important, since although the 
larger collecting area of XMM-Newton, allows for more 
complex spectral analysis, the higher angular resolution 
of Chandra ensures that a source spectrum is not contam- 
inated by diffuse emission and flux from nearby sources. 
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2. OBSERVATIONS AND DATA REDUCTION 

The names of the thirty host galaxies and the details 
of the Chandra ACIS observations are tabulated in Ta- 
ble 1. Distances to the galaxies were obtained from 
ISwartz et al\ (|2004D and references therein This sample 
of gal axies is a subset of those analyzed by ISwartz et all 
( 2004) . No particular criterion was imposed on the selec- 
tion, since our motivation is limited to obtaining enough 
sources and not to study dependency on galaxy type. 

The data reduction and analysis were done using 
CIA03.2 and HEASOFT6.0.2. Using the CIAO source 
detection tool wavdetect, X-ray point sources were ex- 
tracted from the level 2 event list. It was found that at 
least 60 counts are required to fit the spectral data with a 
two parameter model and hence only those sources with 
net counts > 60 were chosen for the spectral analysis. 
Choosing a lower threshold of 50 counts resulted in a 
large number of sources for which spectral parameters 
could not be constrained. To avoid photon pile up ef- 
fects, a conservative threshold of the count rate being 
> 0.05 counts/s was imposed which led to the rejec- 
tion of fifteen sources which have been listed in Table 
3 of the appendix. For some sources, typically near the 
nucleus, it was difficult to find nearby source free back- 
ground regions and hence these sources were also not 
included in the analysis. Sources embedded in exces- 
sive diffuse emission (i.e. when the background flux was 
larger than 2 counts/arcsec 2 ) were also rejected. Typi- 
cally, this amounted to considering only those sources for 
which the estimated background counts were less than 
20. For each data set, observation-specific bad pixel lists 
were set in the ardlib parameter file. Using a combina- 
tion of CIAO tools and calibration data, the source and 
background spectra were extracted. 

These selection criteria makes the sample incomplete 
both in low and high luminosity ends. Thus the results 
obtained should not be used for quantitative analysis of 
the luminosity functions. The motivation here is evaluate 
dependency on spectral models and to identify sources 
which have high intrinsic luminosity. Thus, care has been 
taken to avoid possibly contaminated data, even if such 
criteria result in a loss of sources. 

Spectral analysis was done using XSPEC version 12.2 
and the data was fitted in the energy range 0.3-8.0 keV. 
All sources were fitted with two spectral models, the ab- 
sorbed power-law and an absorbed disk blackbody. Ab- 
sorption was taken into account using the XSPEC model 
phabs. Since the number of counts in each spectrum was 
typically low, the C-statistics was used for the analysis. 
Technically, the C-statistics is not appropriate for high 
counts and/or for background subtracted data. However, 
it was ascertained that the model parameters obtained ei- 
ther by C-statistics or \ 2 statistics, were consistent with 
each other for high count rate sources. This is reassuring, 
because if the results depended on the statistics used, it 
would be imperative to use the correct statistics for high 
count sources (which need not be x 2 ) taking into account 
the correct (possibly non-Gaussian) error profiles. 

An important problem, when fitting low count data 
with a two parameter (plus normalization) model, is 
the possibility of many local minima in the discerning 
statistic (in this case C-statistic) space. Hence, we take 
a cautious approach and do not fit the data using the 




10 3« I I I I I 

10 36 10 3T 10 3a 10 33 10 40 10 41 10 4a 10 43 

Luminsoity L Dbb (ergs/sec) 



Fig. 1. — The intrinsic luminosity (in 0.3-8.0 keV energy range) 
estimated by fitting the power-law model versus the bolometric 
intrinsic luminosity estimated using the disk black body model. 



XSPEC minimization routine. Instead we compute the 
C-statistic for a range of parameter values (using the 
XSPEC command steppar) and find the global minimum. 
Such a technique is numerically expensive, but it ensures 
that the global minimum has been found and the correct 
errors are obtained for the best fit parameters. 

3. RESULTS 

The 365 sources considered in this analysis were classi- 
fied into three categories, depending on whether the data 
was better fit by the disk black body model (23 sources) 
or the power-law one (67 sources) or both (275 sources). 
The criterion chosen to determine a better fit to the data 
was that C-statistic difference between the models should 
be larger than 2.7. If the difference was less, than both 
model fits were considered to equally represent the data. 
Although such a criterion is ad hoc (considering the un- 
certainties in the actual error statistics) and count rate 
dependent, it does serve as an qualitative guideline to 
differentiate between those systems which can be repre- 
sented by a power-law emission and/or a black body one. 
The spectral parameters for all sources for the power-law 
and disk blackbody models are tabulated in Tables 4 and 
5 of the appendix. 

For ULX, the physically relevant parameter is the in- 
trinsic bolometric luminosity which should be used to de- 
fine and identify them. However, given the limitations of 
an instrument's energy sensitivity range, the bolometric 
luminosity is spectral model dependent. For a power-law 
the bolometric luminosity cannot be estimated and only 
a lower limit can be obtained using the observed energy 
range. Since our motivation is to show how the bolomet- 
ric luminosity is affected by the use of different spectral 
models, we have plotted in Figure 1, the bolometric lu- 
minosity for the disk black body model versus the lower 



4 Senorita Devi et al. 



TABLE 1 
Sample Galaxy properties 



Galaxy 


Distance (Mpc) 


ObsID 




N(> 60rfs) 




2.6 


969 


13.98 


13 


IN VjUUUijO 


9 7 


2058 


46.16 


7 


in uVjUoyi 


1 n o 


7Q4 


on 


14 


IN \J V_ _L Ziy _1_ 


8.9 


795 


39 16 


14 


IN uvj -LO-LU 


17.0 


2022 


29.85 


9 


in uvj ioyy 


18.3 


319 


55.94 


36 


NGC1569 


2.2 


782 


96.75 


16 


IN V_xv^Zi 1 ±L/0 


3 1 


2014 


35.59 


4 




3.9 


361 


33.25 


5 


in vju ou i y 


15.6 


2038 


26 57 


5 


in Vj v_ oo i y 


11.1 


1587 


31 52 


7 




14.1 


2025 


59.36 


15 


NGC3628 


10.0 


2039 


57.96 


XI 


NOC4125 


24.2 


2071 


64.23 


8 


IN VjvjiOUU 


20.9 


2015 


40.42 


9 


NGC4374 


17.4 


803 


28.47 


1 


NOC4449 


3.7 


2031 


26.59 


12 


NnP44s t i /on 


7 8 


1579 


19 52 


9 




10.3 


2027 


10.70 




NGC4579 


21.0 


807 


33.90 


3 


NGC4594 


9.6 


1586 


18.51 


18 


NGC4631 


7.6 


797 


59.21 


12 


NGC4649 


16.6 


785 


36.87 


23 


NGC4697 


11.8 


784 


39.25 


19 


NGC5055 


9.2 


2197 


27.99 


16 


NGC5128 


1.0 


962 


36.50 


22 


NGC5194/5 


8.1 


1622 


26.80 


18 


NGC5457 


7.0 


2065 


9.63 


4 


NGC5775 


26.7 


2940 


58.21 


15 


NGC6946 


5.5 


1043 


58.28 


16 



Note. — (T exp )the exposure time in ks; (N) the number of point sources with total counts from the source> 60 as detected by 
wavdetect with fluxscale= 1 



limit to the luminosity using the power-law model. The 
figure represents only those sources which can be rep- 
resented by both models. The figure shows that while 
for most sources the difference in luminosities is not sub- 
stantial, there are sources with estimated luminosities 
> 10 39 erg/s, where the disk black body luminosity es- 
timation is significantly smaller than the power-law one. 
This happens for sources for which the spectral index is 
large. 

To show the dependence of the luminosity function on 
the fitting model, we compute the cumulative luminos- 
ity function in two ways. For the disk black body cu- 
mulative luminosity function (DBCLF), the luminosity 
obtained form fitting a black body is used, except for 
those sources that are fitted better with a power-law 
for which the power-law model estimated luminosity is 
considered. Similarly, for the power-law cumulative lu- 
minosity function, (PLCLF) the luminosity corresponds 
to the power-law fit, except for those sources which are 
better represented by a disk black body spectrum. In 
Figure 2, the solid line represents the DBCLF while the 
PLCLF is plotted as a dotted one. There are less num- 
ber of sources with L > 10 40 ergs/s for disk black body 
preferred representation. Moreover, there is a significant 
difference in the slope of the two luminosity function and 
the presence of a faint break at L ~ 10 40 ergs/s for the 
PLCLF, is not evident in the DBCLF. Of particular in- 
terest is the number of sources whose minimum luminos- 
ity (i.e. the minimum of the two lower limits obtained by 
fitting the two spectral models) exceeds a certain value. 



In Figure 2, the dashed line represents such a minimum 
cumulative luminosity function, which reveals that there 
are eight sources with minimum luminosity greater 10 40 
ergs/s (ELX) and ~ 80 sources with minimum luminosity 
greater 10 39 ergs/s (ULX). 

Figures 3 (a) shows the variation of the luminosity ver- 
sus the disk black body temperature, While most of the 
sources have an inner disk temperature ~ 1 keV, as is 
evident from the distribution (Figure 3 b), there is a pop- 
ulation of high luminosity source with low (~ 0.1 keV) 
temperature. Although the number of sources is low, 
there seems to be some evidence, that ELX (i.e. sources 
with luminosities > 10 40 ergs/s ) can be divided in two 
groups, a "super-soft" group with temperature less than 
0.2 keV and an harder group with temperature ~ 2 keV. 
Figure 3 may be co mpared with the results obtained by 
IWinter et all (|2006f ) using XMM-Newton data for a dif- 
ferent sample and selection criteria. They also find that 
sources with luminosities > 5 x 10 39 ergs/s have have a 
similar bimodal distribution in temperature as shown in 
Figure 3. This supports the hypothesis that ELX can 
be divided into two groups and this is not an artifact of 
sample selection bias. 

Figure 4 (a) shows the variation of luminosity with 
power-law photon index for those sources which can be 
fitted by a power-law model. There is a clear correla- 
tion between the two. This correlation does not seem 
to be due to overestimation of column density, since no 
such correlation is seen in the luminosity versus Nh 
plot (Figure 5 b). Similar to the analysis using the 
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Fig. 2. — The cumulative luminosity function using the disk 
black body model (solid line), the power-law model (dotted line) 
and the minimum cumulative luminosity function (dashed line). 



Fig. 4. — (a) The luminosity versus power-law index and (b) 
the distribution of the power-law index for sources whose spectra 
can be modeled as a power-law emission. The triangles represent 
sources which are better fitted by the power-law model as compared 
to the disk black body one. Two sources which were identified as 
background AGN (see text) are not included in this plot. 
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Fig. 3. — (a) The luminosity versus inner disk temperature and 
(b) the distribution of the inner disk temperature for sources whose 
spectra can be modeled as disk black body emission. The triangles 
represent sources which are better fitted by the disk black body 
model as compared to the power-law one. The two solid lines 
represent the expected luminosity versus maximum temperature 
relations for accretion disks radiating at one and one-tenth of the 
the Eddington Luminosity. Two sources which were identified as 
foreground stars (see text) are not included in this plot. 



disk black body model, there is a group of "super-soft" 
sources (i.e. photon spectral index > 3) which are also 
highly luminous (L > 10 40 ergs/s). The column den- 
sity versus luminosity plots for both the power-law and 
disk black body models (Figure 5), reveal an absence 
of correlation, which is indicative that there may not 
be a bias in the analysis, i.e. the luminosities are not 
being over-estimated because of a Nh overestimation. 
In this analysis, there are eight sources which have an 
apparent minimum luminosity greater than 10 40 ergs/s. 
However, two of these sources (NGC 5055, R.A: 13 15 
30.18, Dec: +42 03 13.5 and NGC 4594, R.A: 12 39 
45.22, Dec: -11 38 49.8) are foreground stars based 
on the optical images of the galaxies. Optical spec- 
troscopy of a source in NGC 5775 ( R.A: 14 53 55.8, 
Dec: +3 33 28.02) reveals that it is a background AGN 
(|Gutierrez fc L6pez-Corredoiral [2005). while a source in 
NGC 139 9 ( R.A: 14 53 5 5.8, Dec: +3 33 28.02) is a 
BLAGN (|Green et al.ll2004D . The spectral properties of 
the other four sources, which we call Extremely Lumi- 
nous sources (ELX) are tabulated in Table 2. The NGC 
0628 source reported in Table 2, is a different source than 
the well studied ULX, CXOU J013651. 1+ 154547. The 
lumin osity of this highly variable source (|Krauss et al.l 
2005) is ~ 10 39 ergs/s and its spectral properties are 
listed in Table 4 . The source in NGC 6946 is a well known 
variable source (|Liu fc M irabel 2 005T) and has been called 
X7 with L ~ lb 39 - 22 e rgs/secdLiraet"al| l20QOh. IXO 85 
dColbert fc Ptakl 120021) and source no. 56 (jHolt et"aL1 
2003). Although, the best fit parameters for C-statistics 
are shown, these results have been checked using % 2 
statistics and by C-statistics fit for unbinned and back- 
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Fig. 5. — The column density versus luminosity for sources fitted 
with (a) the disk black body model and (b) the power-law model. 



ground not subtracted data and it was found that the 
parameters values are consistent within errors and the 
errors on the estimated luminosities do not vary by more 
than a factor two. In general, these sources are better 
represented by disk black body emission than a power- 
law model, except for the source in NGC 6946, which 
however requires an exceptionally large power-law pho- 
ton spectral index, (T > 5). The spectral properties of 
these bright sources suggest that they may be divided 
into two groups. The first group of four sources (Table 
2), are represented by low inner disk temperatures (< 0.5 
keV), and hence may be called "supersoft" sources. In 
contrast the second group of three sources, have harder 
spectra with inner disk temperatures > 1 keV or with 
power-law photon index (T ~ 2) and hence may be called 
hard sources. For this group, the spectra are marginally 
fitted better with a disk blackbody emission, although 
considering the uncertainties in the spectral fitting, a 
power-law representation may also be acceptable. 

4. SUMMARY AND DISCUSSION 

Chandra observations of thirty galaxies were analyzed 
and the spectra of their points sources were fitted us- 
ing both a power-law and a disk black body emission 
model. Only those sources were chosen, which were 
bright enough to allow a meaningful spectral fit, but 
whose data was not contaminated by excessive diffuse 
emission and/or effected by pile-up. It was found that 
the shape of the luminosity function especially at the 
high luminosity end, depends on the choice of the spec- 
tral model. 

In accordance with earlier results, a large number of 
the sources (~ 80 ) have a luminosity which exceeds 10 39 



ergs/s and hence satisfy the standard definition of be- 
ing Ultra Luminous X-ray sources (ULX) and do not 
seem to have any spectral distinction when compared 
with sources having lower luminosity. In this sample of 
365 sources, there are four source which we refer to as 
Extremely luminous X-ray sources (ELX) since their lu- 
minosities were estimated to exceed 10 4 ° ergs/s. These 
sources are in general better described by disk black- 
body emission and can be distinctively grouped into two 
classes. This is consistent with the results of an indepen- 
dent analysis using XMM-Newton data ([Winter et al.l 
2006). The members of the first class have soft spec- 
tra with best fit inner disk temperature < 0.5 keV, while 
for the other class the spectra is harder with inner disk 
temperature > 1.3 keV. 

If disk black body emission is indeed the correct radia- 
tive process for the supersoft class then the inner disk 
temperature should correspond to the maximum color 
temperature of a disk, which can be estimated to be 



T, 



col 



0.3 keV L, 



-1/4 



L.I LEdd 



-1/2 



(1) 



where L40 is the luminosity in 10 40 ergs/s, / is the color 
factor and LEdd is the Eddington luminosity. In Figure 
3 (a), the two solid lines represent this relationship for 
L/ Lsdd = 1 and 0.1. Thus within the uncertainties, the 
supersoft sources are compatible with having pure black 
body disk emission, and have L ~ 0.5LEdd- 

ELX which are members of the hard class have inner 
disk temperatures which are higher than that expected 
from a Eddington limited black body accretion disks. 
Hence, for these source the radiative mechanism is prob- 
ably inverse Comptonization of soft photons. Detailed 
spectral analysis, which included XMM observations, of 
the bright X-ray source in M82 X-l has revealed that its 
spectr um is better fitted by a saturated Comptonization 
model (jAgrawal fc M isra 2006) , which is also the case for 
the the bright X-ray source in Holmberg IX. Holmberg IX 
is not part of the sample studied here and M82 X-l has 
been excluded because of pile-up effects and excess dif- 
fuse emission. Thus, these sources, with estimated lumi- 
nosities ~ 10 41 erg/s, could also be members of the hard 
class of ELX. Thus it seems that like the the hard state 
of standard black hole binaries, the hard class ELX also 
have spectra which is due to thermal Comptonization, 
however unlike black hole binaries, in ELX the Comp- 
tonization seems to be saturated. Thus, it is tempting 
to draw by analogy, that the two spectral classes of ELX 
are actually two spectral states of the same kind of ob- 
ject. This can be verified if spectral transition between 
the two classes is observed. 

With the identification of these ELX and other sources 
from the literature, it is now possible to undertake a more 
extensive study of their properties. Temporal variability 
of these sources will shed more light on the nature of 
these enigmatic sources. 

The authors thank the referee for useful comments and 
suggestions which have significantly improved the paper. 
ASD thanks CSIR and IUCAA for support. 
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Spectral properties of sources with minimum luminosity > 10 40 ergs/s 
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log L Pow 


C _ Str£ltr PoW 
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NGC6946 
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20 35 0.13 


+15 47 45.01 
+-60 9 7.97 


200 
1936 


7 

66 


09+ ' 04 
u.uy_ 03 

n o 9 +0.02 
u - oz -0.03 


41.65l?;S 
39.461°°* 


0.7 
162.2 


5-i4±g;|* 


42.18+?'^ 


8.6 
134.8 


NGC4579 
NGC5775 


12 37 40.30 
14 53 58.90 


+11 47 27.48 
+ 3 32 16.78 


1696 
1358 


66 
60 


i oc+0.10 
1 - oo -0.09 

1.97±°;g 


40.56tg;g 

4i.n±8:S 


65.0 
55.8 


i 7 r+0.11 

l -'°-o.u 
1 86+ ' 22 


40.371°;°° 
40.95l°; ° 


74.5 
60.1 



Note. — Host Galaxy name; Right Ascension; Declination; Total counts; Number of energy bins after rebinning; Best fit inner disk 
temperature; Bolometric luminosity estimate using disk black body model; C-statistic for disk black body model; Photon index T for the 
power-law model; Luminosity estimate (0.3-8.0 keV) using the power-law model; C-statistic for power-law model; 
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TABLE 3 

List of sources with count rate greater than 0.05 counts/s 



Galaxy 


R.A. 


Dec. 


Count rate s 1 


JNGC0253 


r\ at on at 

47 32.97 


or" 1 T AO O O 

-25 17 48.80 


O O O A F" 1 n A 

0.0845164 


NGC0253 


47 22.59 


-25 20 50.87 


0.0657110 


NGC0253 


/~v AT IT f P" 

47 17.55 


or" "io 11 10 

-25 18 11.18 


O OF'/^OOOT 

0.0569807 


in LjCiooy 


4 ol 10. OO 


+04 4y OU.lo 


U.UOO 1 4oo 


NGC2403 


7 36 55.61 


+65 35 40.85 


0.0746877 


NGC2403 


7 36 25.53 


+65 35 40.02 


0.1521207 


NGC3628 


11 20 15.75 


+13 35 13.70 


0.0530390 


NGC4374 


12 25 11.92 


+12 51 53.53 


0.0563400 


NGC4449 


12 28 17.83 


+44 6 33.86 


0.0525013 


NGC4485 


12 30 43.26 


+41 38 18.36 


0.0509909 


NGC4485 


12 30 30.56 


+41 41 42.33 


0.0758634 


NGC4559 


12 35 58.56 


+27 57 41.91 


0.1233310 


NGC4559 


12 35 51.71 


+27 56 4.05 


0.1984267 


NGC4631 


12 41 55.56 


+32 32 16.90 


0.0554584 


NGC6946 


20 35 0.74 


+60 11 30.74 


0.1448118 



Note. — The spectra of these sources would be affected by pile-up and hence have not been included in the sample 



TABLE 4 

Spectral Properties of point sources fitted with the Power-Law model 



Galaxy 


R.A. 


Dccl. 


n H (10 22 cm- 2 ) 


r 


log(L) crgs/s 


Cstat 


d. 0. f. 


NGC0253 


47 43.07 


-25 15 29.28 


15+ - 53 
u - io -0.l5 


o cn+2.31 
z *° — 1.21 


37.70± ; 29 


9.43 


4 


NGC0253 


47 42.80 


-25 15 2.02 


0-77t°oj 9 


1.86+g- f 7 


38.26 + °- 2 2 


3.36 


4 


NGC0253 


47 35.25 


-25 15 11.53 


o.48± ;?! 


2 1Q+° :22 


38.651°;°! 


38.94 


19 


NGC0253 


47 34.28 


-25 17 3.32 


5 66+ 5 - 98 

0.00_ 5 gg 


t- on+4.62 
J.OO_ 4 .g 5 


40.82l 3 2 ; 28 


3.21 


2 


NGC0253 


47 34.00 


-25 16 36.51 


1 04+ 011 


z - la -0.22 


on n1 +0.08 


10.93 


27 


NGC0253 


47 33.55 


-25 18 16.51 


oo+° 18 


1 K4+0-99 
i - o4 -0.44 


07 00+0.15 
ot.oo_ Q 10 


2.37 


2 


NGC0253 


47 32.05 


-25 17 21.43 


o -, 7 +2.20 
-2.46 


J-- 9 ' -1.10 


38-6il8:g 


1.63 


3 


NGC0253 


47 30.98 


-25 18 26.23 




1.751111 


38.24+.°;^ 


1.76 


3 


NGC0253 


47 28.01 


-25 18 20.21 


0-"3_ ; 92 


2.08ti;™ 


38.02±J;§| 


2.13 


1 


NGC0253 


47 25.20 


-25 19 45.22 


00+° 15 

u - uu -o.oo 


-, 01+O.66 
1 - O1 -0.22 


37 - 96 -0.09 


5.10 


4 


NGC0253 


47 18.50 


-25 19 13.94 


00+° 15 
u - uu -o.oo 


-, O-.+0.44 
1 - O1 -0.22 


07 q 7 +0.07 
— 0.05 


3.94 


5 


NGC0253 


47 40.66 


-25 14 11.71 


n 69+°' 72 


2 41+ 1 - 98 
^• 4i -1.32 


37.96± ; 27 


1.07 


3 


NGC0253 


47 17.65 


-25 18 26.45 


n+ - 49 
u - il -o.n 


1 01 +1.54 


38.121 ,; 27 


3.13 


3 


NGC0628 


1 36 51.06 


+15 45 46.86 


0.03±o.05 




39 22+ 04 


50.58 


36 


NGC0628 


1 36 47.45 


+15 47 45.01 




9.56_ 2 42 


42.18t°-j6 


8.56 


4 



Note. — Host galaxy name; Right Ascension; Declination; nu, equivalent hydrogen column density; T, photon power-law index; L, 
X-ray luminosity in the energy range: 0.3-8.0 keV; C-statistics; degree of freedom. The complete version of this table is in the electronic 
edition of the Journal. The printed edition contains only a sample. 
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TABLE 5 

Spectral Properties of point sources fitted with the disk black body model 



Galaxy 




R.A. 


Decl. 


NGC0253 





47 


43.07 


-25 15 29.28 


NGC0253 





47 


42.80 


-25 15 2.02 


NGC0253 





47 


35.25 


-25 15 11.53 


NGC0253 





47 


34.28 


-25 17 3.32 


NGC0253 





47 


34.00 


-25 16 36.51 


NGC0253 





47 


33.55 


-25 18 16.51 


NGC0253 





47 


32.05 


-25 17 21.43 


NGC0253 





47 


30.98 


-25 18 26.23 


NGC0253 





47 


28.01 


-25 18 20.21 


NGC0253 





47 


25.20 


-25 19 45.22 


NGC0253 





47 


18.50 


-25 19 13.94 


NGC0253 





47 


40.66 


-25 14 11.71 


NGC0253 





47 


17.65 


-25 18 26.45 


NGC0628 


1 


36 


51.06 


+15 45 46.86 


NGC0628 


1 


36 


47.45 


+15 47 45.01 



n H (10 2 



kT ln (keV) 



log(L) ergs/s 



d. o. f. 



00+ ' 32 
u - uu -0.00 

44+ - 36 
u -^-0.32 

24+ - 07 
u - z4 -0.07 
R +5.36 
-3.07 
,,+0.12 
-0.10 
-1+0.08 
-0.00 
hi. 82 
-1.62 
,+0.96 
-0.86 

62+ - 77 
u - D ^-0.56 
-,+0.06 
-0.00 
-1+0. 05 
-0.00 

42+ - 47 

-0.38 

o.08tjj;g? 

00+ 00 
-+0.58 

0.39 



3.08Z 
0.64^ 

o.ooi* 

2.04j 
0.96 4 



o.oo: 

0.00 4 



0.751 



0.59 
1.56 



-0.25 
-1.95 
-0.52 
+0.21 



„+0.19 
-0.16 
-+0.69 
3 -0.29 
1-4.71 

-o.so 

B.00 
-0.91 
1+8.80 
-0.55 

1 20+ 086 
± - zu -0.33 
3+0.73 
-0.41 

qn+ 2 ' 76 

u - au -0.40 

1 24+ 8 ' 76 
1 -^ 4 -0.65 

1+0.08 
'-0.11 
1+0.04 
-0.03 



i.34: 

0.! 
2.07^ 
2.00^ 
1.20" 1 



1.59"* 



O.i 

0.09^ 



37.69t°;g 7 
38.40±°;j» 

,+0.04 
L -0.03 

38.69± 2 ;°S 
39.0610^3 

o 7 40+O.2O 

1-0.40 
-0.09 
,+0.73 



38.7V 



38.68^ 



38.411 012 
38.081°;?° 

00 1 -21 

00.1Z_Q Q 9 

38.22_ ; og 

37 q.4+0-38 

01 - y -0.18 

OO QQ + 1.04 

— 13 

39.42+H 2 
41.651?$ 



9.41 
3.95 

35.79 
3.37 
8.60 
3.85 
1.37 
1.84 
2.51 
7.71 
6.36 
1.64 
2.85 

119.44 
0.66 



4 
4 

19 
2 

27 
2 
3 
3 
1 
4 
5 
3 
3 

36 
4 



Note. — Host galaxy name; Right Ascension; Declination; uh , equivalent hydrogen column density; kTi n , inner disk temperature; L, 
Bolomctric X-ray luminosity; C-statistics; degree of freedom. The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 



